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We present here a nonlinear sigma-model (NLcrM) study of a spin-1 antiferromagnetic Heisenberg 
chain in an external commensurate staggered magnetic field. We find, already at the mean-field 
level, excellent agreement with recent and very accurate Density Matrix Renormalization Group 
(DMRG) studies, and that up to the highest values of the field for which a comparison is possible, 
for the staggered magnetization and the transverse spin gap. Qualitative but not quantitative 
agreement is found between the NLcrM predictions for the longitudinal spin gap and the DMRG 
results. The origin of the discrepancies is traced and discussed. Our results allow for extensions 
to higher-spin chains that have not yet been studied numerically, and the predictions for a spin-2 
chain are presented and discussed. Comparison is also made with previous theoretical approaches 
that led instead to predictions in disagreement with the DMRG results. 

PACS numbers: 75.10.Jm, 75.0.Cr, 75.40.Cx 



One and quasi-one-dimensional magnets (spin chains 
and ladders) have become the object of intense analyt- 
ical, numerical and experimental studies since Haldane 
put forward |JjJ his by now famous "conjecture" according 
to which half-odd-integer spin chains should be critical 
with algebraically decaying correlations, while integer- 
spin chains should be gapped with a disordered ground 
state HII- Haldane's conjecture has received strong ex- 
perimental support from neutron scattering experiments 
on quasi-one-dimensional spin-1 materials such as NENP 
and Y2BaNiC>5 0]. It is by now well established that 
pure one-dimensional Haldane systems (i.e. integer-spin 
chains) have a disordered ground state with a gap to a 
degenerate triplet (for spin-1 systems) of magnon exci- 
tations. There is also general consensus on the value of: 
Ao = 0.41048(2)J (with J the exchange constant) for 
the gap in spin-1 chains, that has been obtained by both 
Density Matrix Renormalization Group (DMRG) p] and 
finite-size exact diagonalization [E§ methods. 

The effects on Haldane systems of external magnetic 
fields have also been the object of intense studies, again 
both experimental, numerical and analytical. A uni- 
form field induces a Zeeman splitting of the degener- 
ate magnon triplet, with the Haldane (gapped) phase 
remaining stable, with zero magnetization, up to a lower 
critical field H Cl = Ao, when one of the magnon branches 
becomes degenerate with the ground state and a Bose 
condensation of magnons takes place j7Uq]. At higher 
fields the system enters into a gapless phase and magne- 
tizes, with the magnetization reaching saturation at an 
upper critical field H C2 = AS J . 

In view of the underlying, short-range antiferromag- 
netic (AFM) ordering that is present in the Haldane 
phase, the study of the effects of a staggered magnetic 



field appears to be even more interesting. Of course 
static, staggered fields cannot be manufactured from 
the outside. However, recently a class of quasi-one- 
dimensional compounds has been investigated that can 
be described as spin-1 chains acted upon by an effective 
internal staggered field. Such materials have Y2BaNi05 
as the reference compound and have the general for- 
mula R2BaNi05 where R (^ Y) is one of the magnetic 
rare-earth ions. The reference compound is found to 
be highly one-dimensional, hence a good realization of 
a Haldane-gap system (remember that the Ni 2+ ions 
have spin S = 1) ||. The magnetic R 3+ ions are po- 
sitioned between neighboring Ni chains and weakly cou- 
pled with them. Moreover, they order antiferromagneti- 
cally below a certain Neel temperatute T/v (typically: 16 
K< T N < 80 K ||). This has the effect of imposing an 
effective staggered (and commensurate) field on the Ni 
chains. The intensity of the field can be indirectly con- 
trolled by varying the temperature below T/\r. The stag- 
gered field lifts partially the degeneracy of the magnon 
triplet, leading to different spin gaps in the longitudinal 
(i.e. parallel to the field) and transverse channels. Neu- 
tron scattering experiments on samples of Nd2BaNi05 
and Pr2BaNi05 |l(| show an increase of the Haldane gap 
as a function of the staggered field. Experimental results 
have also been obtained for the staggered magnetization 

& 

Recently, an extensive DMRG study of an S = 1 
Heisenberg chain in a commensurate staggered field has 
appeared |l^| . The authors in Ref . |l^] have obtained ac- 
curate results for the staggered magnetization curve, the 
longitudinal and transverse gaps and the static correla- 
tion functions (both longitudinal and transverse). They 
found however a strong disagreement in the high-field 
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regime between their results and previous theoretical ap- 
proaches |l3|,[l4j whose starting point was a mapping of 
the chain onto a nonlinear sigma-model, and this has led 
them to openly question the entire validity of the NLctM 
approach, at least in the high-field regime. 

This strong criticism has partly motivated our study 
of the same model as in Ref. . We will actually show 
that an accurate treatment of the NLcrM does indeed 
lead to an excellent agreement between the analytical 
and the DMRG results both for the magnetization and 
the transverse gap for all values of the staggered field. 
Some discrepancies are instead present for the longitu- 
dinal gap between our results and the data of Ref. [03 , 
but our identification of the latter relies on an approxi- 
mation that, as we shall argue in the concluding part of 
this Letter, may be questionable in the presence of an 
external (staggered) field. After presenting the deriva- 
tion of our results, we will also discuss what are the rea- 
sons for the disagreement between previous theoretical 
approaches and the DMRG results. 

We start from the following Hamiltonian for a Heisen- 
berg chain coupled to a staggered magnetic field: 



n 



H = 5^[JS i .S<+i+H..(-l) i S i ] 



(1) 



where Sf = S(S + 1) (we set h = 1 from now on), the 
staggered magnetic field H s includes Bohr magneton and 
gyromagnetic factor, and J > is the exchange constant. 
In this Letter we will consider the case of integer spin 
S. Under this assumption, the Haldane mapping of the 
Heisenberg chain onto a (1+1) nonlinear (T-model can 
be readily obtained |l4j]. The topological term is absent 
since the spin S is integer, while the staggered magnetic 
field couples linearly to the NLcrM field n. The euclidean 
NLcrM Lagrangian is thus given by E3| : 



The euclidean action Se is quadratic in the field n 
which can be integrated out exactly. To this end, we 
introduce the notation x= (x,t) and write 



S E = J dxdx' JT(x,x')n(x) . n(x') 
- SH S ■ Jdxn(x)+i J dxA(x) , 



(4) 



where the kernel K is given by 
1 



K(x,x') 



2gc 



c 2 dl + d 2 T + 2igc\{x)]8{x-x') . (5) 



The square in equation (|J) can be completed by perform- 
ing the linear shift n(x) — > n'(x) = n(x) — n(x), where 



S 



H 



(6) 



and K~ x is the inverse of the kernel (||). After the in- 
tegration over the shifted field we obtain the effective 
action for the field A 



5(A) = ^TrlnX + i J rfxA(x) 

\ J dxdx' J fir- 1 (x,x') 
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(7) 



The partition function is now given by Z = f[^-]e~ s ^' 
and can be calculated by a saddle-point approximation 
for the integral over A. We look for a static and constant 
saddle-point solution: the kernel K depends in this case 
only on the relative coordinates K(x,x') = K(x — x'). 
reads then: 



The equation 



-A'- 1 (x = 0) = l-^^ 
2 y ' 4 



J dxK- l {x) 



(8) 



C E = — 



1 
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On 
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S 1 -*' Su ' H * 



(2) 



where the NLcrM field n is a three-component unit vec- 
tor pointing in the direction of the local staggered mag- 
netization. The bare coupling constant g and spin-wave 
velocity c are related to the parameters of the Heisen- 
berg Hamiltonian by g = 2/S and c = 2JSa, with a 
the lattice spacing. The constraint n 2 = 1 can be taken 
into account in the path-integral expression for the par- 
tition function by writing Z = f [2?n] [Ip] e _s - E , where 

S E = dx Jq drC(x, t) [L = Na being the total length 
of the chain) and 

C(x, t) = Ce(x, t) — iX(x, r)[n(x, r) 2 — 1] . (3) 

The functional integration over the field X(x,t) imple- 
ments the constraint n 2 = 1. 



It is convenient to work now in Fourier space, which is 
defined according to 



n(q,n) 



L f f3 

dx / dTe- t{qx -^n{x,T) 



(9) 



(the frequencies 0„ = 2-nnj ft are Matsubara Bose fre- 
quencies), and similarly for the field A and kernel K. 
The saddle-point equation becomes: 



" £*-%n) = l-^ 



2(pL) 



K- l {q,n) 



2gc 



f2 2 + c 2 c/ 2 - i\2gc 



K- 1 (q = 0,n = 0) (10) 
(11) 



We now make the assumption —i\2gc = c 2 /£ 2 , with £ 2 
real and positive. Under this assumption, that will be 
consistently verified shortly below, the sum over frequen- 
cies in Eq. ( |To| ) can be performed by standard techniques. 
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By taking the thermodynamic and zero-temperature 
limit (L, (3 — > oo) we finally obtain the following equation 
for £ 



2?r 



AC + v/i + (A|p 



= 1 - 



„9 S 



(12) 



where we have introduced the ultraviolet cutoff A to reg- 
ularize the momentum integration. The staggered mag- 
netization (in units of S) will be determined instead by 
the equation for the linear shift, which now reads 



K~ x {q = 0,n = 0) = 



gse 



(13) 



In order to get rid of the ultraviolet cutoff A we observe 
that from Eq. ( |Tl| ) it follows that, when H s = 0, the 
Green functions for the field n are given by 



G a ,p{q,n) 



h a (q,n)h f3 (-q, - 

gc 



-n)) 



8 a , 



c 2 q 2 



2 /e 



(14) 



with £ = sinh(27r/3g)/A, as obtained from Eq. ( p2[ ) for 
H s = 0. From the structure of the Green functions Jl4| ) 
it follows that the excitations are gapped, with the value 

A = A(H S = 0) = . - ^ C (15) 

for the zero-field Haldane gap |jj Aq. The Haldane gap 
is the only free parameter in our theory, that we will as- 
sume from previous numerical estimates |5],[3| . In this way 
we fix the cutoff. Eqs. (|lj) and (|l^) (with the addition 
of Eq. (|l5|)) provide our staggered-magnetization curve. 
Eq. ([[2]) can be solved analytically for large or small val- 
ues H s , while in general it has to be solved numerically 
(albeit with no effort). 
For small EL we obtain 



m. 



gSc 
A 2 



H, 



1 - 



4tt/3 c 2 gS 2 
tanh(27r/3g) A 4 



(16) 



The zero-field staggered susceptibility is then given by 



X.(0) 



dm, 



H s =0 



gSc 
Af" 



4JS 

A 2 
^0 



(17) 



For S = 1, by taking A = 0.41048 J from Ref. f§, we 
obtain x s (0) = 23.74/J which agrees fairly well with 
Monte-Carlo result x*(0) = 21(1)/J (l| and DMRG 
calculation Xs(0) = 18.5/J @. For S = 2, we take 
A = 0.0876J from Ref. @ and get % s (0) = 1043/ J, 
which should be checked by new numerical investigations. 

For large H s the staggered magnetization is instead 
given by 



A 



+ 



2H< 



0(A 3 /H^ /2 ) 



(18) 




0.0 0.1 0.2 0.3 0.4 0.5 

FIG. 1. The staggered-magnetization curve for S = 1. Our 
results (solid line) are compared with DMRG data (circles) 
and with the NLctM results (dotted line) of Refs. @ and @, 
respectively. 



where A = A ^ *j00$ ■ It is obvious from Eq. © 
that the magnetization saturates to its maximum value 
of m s = 1 only asymptotically for H s — > 00. This is 
consistent with the fact that it is only in that asymptotic 
limit that a fully polarized Neel state becomes an exact 
eigenstate of the original Heisenberg Hamiltonian (Q) , as 
well as with the fact that our mean-field approach keeps 
track of the NLcM constraint, albeit only on the average. 

In Figure 1 we compare, for S = 1, our staggered- 
magnetization curve with DMRG data |L2l and with the 
NLcrM treatment of Ref. |l4| . The magnetization curve 
was obtained by solving numerically Eq. (|l2|) (with A 
given by Eq. Jl5|)) and by inserting the resulting £ into 
Eq. ([lj) . As already anticipated, the agreement with the 
DMRG data is excellent over the whole range of fields. 

We have also solved Eq. ([l2|) for S = 2, by fix- 
ing again A via Eq. (n5| with the value, specific to 
S = 2, A = 0.0876 J |n|. At present, the staggered- 
magnetization curve for S = 2 has not yet been obtained 
by numerical methods. Given the excellent agreement 
between our curve and numerical data for S = 1, we 
expect that our magnetization curve (shown in Figure 
2) will agree with future numerical investigations also 
in the case S = 2. Indeed, on general grounds, Hal- 
dane's mapping is expected to work even better when 
the value of the spin S is increased. Note that the mag- 
netization increases faster with H s (or, equivalently, the 
staggered susceptibility becomes larger) when the spin S 
is increased. This is consistent with the fact that, when 
S —>■ 00, the system approaches the classical one, which 
is critical at T — and H s = 0. Note also that the value 
of the zero-field Haldane gap roughly sets the scale of the 
staggered-magnetization curve. 
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0.00 0.01 0.02 0.03 0.04 0.05 
FIG. 2. The staggered-magnetization curve for S=2. 



Introducing an external space and time-dependent 
source field J = J(x) and promoting the partition func- 
tion Z to a generating functional Z[3] obtained by re- 
placing in the original path-integral the Euclidean La- 
grangian of Eq. (0) with: 



2ff 



1 fdn 



dn 

Ox 



-n(x)-J(x) , (19) 



the connected Green functions for the staggered field n: 
G^(x,x') = (T r [n Q (x)r/(x')]) - (n"(x)>(/(x')) (20) 

can be obtained by double functional differentiation with 
respect to the source field as: 



G^(x,x') = 



5 2 lnZ[J] 



(21) 



J=SH a 



The saddle-point analysis can be performed in this case 
as well. The saddle-point equation reads now: 

^(x.x) = 1 - I J dydy'K- l (y,x)K- l (x,y>) 

xJ(y)-J(y') . (22) 

This equation will yield in general a space and time- 
dependent saddle point X sp — A sp (x), which will reduce 
however to a constant for J = const. = 5H S when, as it 
can be checked easily, Eq. ( p2| ) reproduces our previous 
Eq. (||). What matters here is the fact that the sad- 
dle point will depend quadratically on the source field. 
From this it follows that the transverse Green functions 
G xx (it, x') = Gj,j,(x, x') do not get any additional contri- 
bution from the dependence of A sp on J (we assume here 
Hs = zH s ) and can be calculated by taking J = SH S 
from the beginning. The transverse Green functions are 
thus still given by the same expression valid for H s = 0, 
Eq. (Jul) where now £ depends on H s via Eq. (OJ). The 



transverse gap will be accordingly Aj_(H s ) 



?(H,) 



, and 



in Figure 3 we compared the transverse gap as calculated 
by this relation with the one obtained by DMRG data. 
The agreement is again excellent. 



< 
< 




FIG. 3. Plots of the transverse (solid curve) and longitu- 
dinal (dashed curve) gap predicted by the NLaM approach, 
compared with the corresponding DMRG data (circles and 
triangles, respectively). 

In the parallel channel, on the other hand, the depen- 
dence of A sp on J yields additional contributions to the 
Green functions, which do not allow us to calculate the 
longitudinal Green functions explicitly. To obtain an esti- 
mate for the longitudinal gap we resorted to the so called 
"single- mode approximation" []l4|| , that is we assumed 
the gap and susceptibilities to be connected by the same 
law both in the transverse and parallel channel. It can 
be readily shown that for an isotropic system, for which 
the magnetization is directed along the external field, 
X± = mf = m s/H s , while x\\ = afff = sjf s - From 
Eq. (JTJ) we thus find X ± = Sg£ 2 /c = Sgc/{A ± ) 2 . If wc 
assume the same relation to hold also in the parallel chan- 
nel, A || can be calculated as a derived quantity from the 
magnetization curve by using A|| = (Sgc/xw) 1 ^ 2 ■ Our 
data for Ay are also shown in Figure 3. The agreement 
with the DMRG data is definitely worse in this case, and 
this may attributed to a shortcoming of the single-mode 
approximation in the longitudinal channel. The complete 
structure of the Green functions along the lines outlined 
above in Eqs. ([l9|) and (21) is being presently under- 



taken, and more extended and complete results will be 
published elsewhere p7[ . 

In conclusion, we come to a brief comparison between 
our approach and that of Refs. jl3| and [ O ]. 

What the two approaches have in common is the start- 
ing point, namely the mapping of the original problem 
onto a NLerM, slightly modified by the external stag- 
gered field. In the abovemntioned references, however, 
the authors resolved to relaxing the NLtrM constraint 
by replacing it with a polynomial self-interaction of the 
n-ficld, keeping terms in the expansion of the interac- 
tion up to eight order. This leads to a (generalized) 
Ginzburg-Landau-type theory with up to six free param- 
eters which the authors took from previous numerical 
and Rcnormalization-Group studies. Their results for 
the staggered magnetization saturate at a finite value of 
the staggered field (see Figure 1), which indicates that 
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softening the NLerM constraint is inappropriate at high [17] E.Ercolessi, G.Morandi, P.Pieri and M.Roncaglia, in 

fields. Also, their results for the gaps deviate badly from preparation. 

the DMRG results. All this has been evidenced in Ref. 

Jl2| and justifies the negative comments of the authors in 

that reference, that however should be addressed more to 

the additional approximations adopted in Refs. ]l3| ] and 

fbij than to the NLerM approach itself. 
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